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The field equations in f(R) gravity derived from the Palatini variational principle and formulated 
in the Einstein conformal frame yield a cosmological term which varies with time. Moreover, they 
break the conservation of the energy-momentum tensor for matter, generating the interaction be¬ 
tween matter and dark energy. Unlike phenomenological models of interacting dark energy, f(R) 
gravity derives such an interaction from a covariant Lagrangian which is a function of a relativis- 
tically invariant quantity (the curvature scalar R). We derive the expressions for the quantities 
describing this interaction in terms of an arbitrary function f(R), and examine how the simplest 
phenomenological models of a variable cosmological constant are related to f(R) gravity. Particu¬ 
larly, we show that Ac 2 = H 2 ( 1 — 2 q) for a flat, homogeneous and isotropic, pressureless universe. 
For the Lagrangian of form R — 1 /R, which is the simplest way of introducing current cosmic ac¬ 
celeration in f(R) gravity, the predicted matter-dark energy interaction rate changes significantly 
in time, and its current value is relatively weak (on the order of 1% of Ho), in agreement with 
astronomical observations. 

PACS numbers: 04.50.+h, 95.36.+X, 98.80.-k 


I. INTRODUCTION 

Einstein’s general relativity is based on the Lagrangian 
that is a linear function of the Riemann curvature scalar. 
The linearity of the gravitational action with respect 
to the curvature results in the Einstein equations of 
the gravitational field that are second-order differential 
equations for the metric tensor ,l|. The Einstein equa¬ 
tions applied to a homogeneous and isotropic spacetime 
lead to the Friedmann equations describing an expand¬ 
ing universe with a positive deceleration. The most 
accepted explanation of the observed cosmic accelera¬ 
tion 0 is that the universe is dominated by dark energy 
(quintessence) 00 The simplest way of introducing 
dark energy into general relativity is to add a cosmo¬ 
logical constant A to the curvature scalar R in the La¬ 
grangian (the AC DM model): 

s g = -^- c j d 4 x^g[ R + 2A}. (1) 

However, it is also possible to modify the curvature 
dependence of the gravitational action to obtain the field 
equations that allow accelerated expansion. A particular 
class of alternative theories of gravity that has recently 
attracted a lot of interest is that of the f{R) gravity mod¬ 
els, in which the gravitational Lagrangian is a nonlinear 
function of R [SJ. It has been shown that current cosmic 
acceleration may originate from the addition of a term 
R _1 (or other negative powers of R) to the Einstein- 
Hilbert Lagrangian R @J. As in general relativity, f{R) 
gravity theories obtain the field equations by varying the 
total action for both the field and matter. 

In the literature, there are two approaches on how to 
perform the variation. We use the metric-affine (Pala¬ 
tini) variational principle, according to which the metric 
and connection are considered as geometrically indepen¬ 
dent quantities, and the action is varied with respect to 


both of them. The other one is the metric (Einstein- 
Hilbert) variational principle, according to which the ac¬ 
tion is varied with respect to the metric tensor g M!/ , and 
the affine connection coefficients are the Christoffel sym¬ 
bols of g^. The field equations in the Palatini formalism 
are second-order differential equations, while for metric 
theories they are fourth-order. Both approaches give the 
same result only if we use the standard Einstein-Hilbert 
action (7]. 

An interaction between ordinary matter and dark en¬ 
ergy has been introduced in the form of a time depen¬ 
dent cosmological constant A (t) flj. A decreasing (de¬ 
caying) A gives a positive rate of the entropy produc¬ 
tion and could explain the observed large entropy of the 
universe. The cosmological constant can depend on the 
cosmic time directly (9j or through other cosmological 
variables such as the scale factor a or the Hubble param¬ 
eter H. Dimensional arg uments lead to A ex a~ 2 jlfilllj . 
A oc H 2 0 , or both llfl . A more general form A ex a m , 
where m = const ]lj|, is predicted by holographic cos¬ 
mology 000 The cosmological term can also in¬ 
volve other quantities such as the deceleration parameter 
or temperature. Ref. 0] lists and reviews phenomeno¬ 
logical models with a variable cosmological term, and 
examines evolution of the scale factor in these models. 

The field equations in f(R) gravity formulated in the 
Einstein conformal frame not only yield an effective cos¬ 
mological term, but also break the conservation of the 
energy-momentum tensor for matter jl9j , generating the 
interaction between matter and dark energy. Therefore, 
f(R) gravity provides a physical explanation for the ex¬ 
change of energy between cosmological term and other 
forms of matter. Unlike phenomenological models of 
interacting dark energy, and like Brans-Dicke cosmolo¬ 
gies with the cosmological term depending on the scalar 
field [20(l . it derives such an interaction from a covari¬ 
ant Lagrangian in which the cosmological constant is a 
function of a relativistically invariant quantity ( R). Fur- 
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thermore, f(R) gravity in the Einstein frame generates a 
variation of the gravitational constant with the cosmolog¬ 
ical time, as in the large number hypothesis of Dirac [2l[ . 

The aim of this paper is to find how this interaction 
depends on the function f{R). In Sec. II we derive the 
equations of field in f{R) gravity using the Palatini vari¬ 
ational principle. In Sec. Ill we apply these equations 
to a homogeneous and isotropic, flat universe filled with 
dust, and obtain the expressions describing the interac¬ 
tion between matter and dark energy. We also exam¬ 
ine how the simplest phenomenological models of a vari¬ 
able cosmological constant are related to f(R) gravity. 
In Sec. IV we give numerical predictions for the case 
f(R) = -R + const. x R _1 , which is the simplest way of in¬ 
troducing current cosmic acceleration in f(R) gravity ($], 
and compare them with other models of interacting dark 
energy. The results are summarized in Sec. V. 


II. THE FIELD EQUATIONS 

The action for an f(R) gravity is given by mu] 

Sj = / d HVHf(R)\ +Sm(g^,ip). ( 2 ) 

Here, \f—gf(R) is a Lagrangian density that depends on 
the curvature scalar R = R Atl/ (r p A cr )g M1/ , S m is the action 
for matter represented symbolically by 1 p and assumed 
to be independent of the connection r p A CT ), and the con¬ 
nection is assumed to be symmetric (no torsion). The 
variables describing an f(R) Lagran gian are said to form 
the Jordan conformal frame ( JCF) |2J|. Tildes indicate 
quantities calculated in this frame, e.g., is the JCF 
metric tensor. 

Variation of the action Sj with respect to g yields 
the field equations 

f'{R)R„u -\f(R)~g^ = kT^, (3) 


does not contain torsion [23] ■ Such a transition is also 
possible for more general theories in which the gravita¬ 
tional Lagrangian depends on the Ricci tensor, but not 
on the Weyl tensor j24|. In the Jordan frame, the connec¬ 
tion is metric incompatible unless f(R) = R. We regard 
the ECF metric tensor as physical, although whether it is 
true or not should be ultimately decided by experiment 
or observation. Since the physical equivalence between 
both frames is an open problem as well |25|, so is the 
answer to the question whether and how the results of 
this paper would change if the JCF was physical. 

In the Einstein frame, the action © becomes the 
standard general-relativistic action of the gravitational 
field interacting with an additional nondynamical scalar 
field [13 El: 


Se = -^~ c j #xy/=g[R-2V{R)\ 
+5 m ([/'(R)]- 1 5 ^,V’), (6) 

where R = R M1/ (r p A 0 .)g MJ ' is the Riemannian curvature 
scalar of the metric g /lv , and V(R) is the effective poten¬ 
tial, 


V{R) = 


Rf'(R)-f(R) 
2 if'(R)} 2 


( 7 ) 


The curvature scalars in both frames are related by 


R = f\R)R , (8) 


which follows from ©• 

Variation of the action m with respect to g^ v gives 
the equation of the gravitational field in the Einstein 
frame mil: 

1 kT 

R, v --Rg^ = -^--V(R)g, v . (9) 

z J \R) 

Eqs. © and © yield an algebraic relation 


where the dynamical energy-momentum tensor of matter 
T pj , is generated by the JCF metric tensor: 

6Sm = h J (4) 

The prime denotes the derivative of the function f(R) 
with respect to R. From variation of Sj with respect to 
the connection Y£ v it follows that this connection is given 
by the Christoffel symbols of the conformally transformed 
metric JB] 


Rf'(R) - 2 f(R) = KTf'(R), (10) 

from which we obtain R as a function of T = for 

a given f(R) [Hj. If T = 0 (vacuum or radiation) then 
R = const, and the solution of the field equation is an 
empty spacetime with a cosmological constant [jj. Sub¬ 
stituting R into © leads to a final relation between the 
geometrical tensors and the energy-momentum tensor, 

R^u - \Rg^u = «r(F)T p! , + A (T)g^, (11) 


g^v — f ( R)g /ji/- (5) 

The metric g^ v defines the Einstein conformal frame 
(ECF), in which the connection is metric compatible [23j. 

A transition from the JCF (modified gravity) to the 
ECF (general relativity) in f(R) gravity is possible if S m 


with a running gravitational coupling n r (T) = 
k[/'(R(T))] _ 1 and a variable cosmological term A (T) = 
-V(R(T)): 


A (R) 


f(R) - Rf'(R) 
2 lf'(R)} 2 


( 12 ) 





3 


Such a relation is, in general, nonlinear and depends on 
the form of the function f(R). Moreover, the time de¬ 
pendence of T makes the quantities n r and A variable, 
although they are not associated with any dynamical 
quintessence scalar field as in 10 . Concluding, Pala¬ 
tini f(R) gravity does not change the structure of Ein¬ 
stein’s general relativity but introduces a varying grav¬ 
itational constant [ 2 l|, as well as a term acting like a 
cosmological constant which varies with time. It modi¬ 
fies the gravitational field inside material objects, while 
the metric in vacuum (such as the Schwarzschild solution 
for spherically symmetric systems) is the same as that in 
general relativity. 


III. THE INTERACTION BETWEEN MATTER 
AND DARK ENERGY 


The Bianchi identity applied to © gives 


T„ 



7W) 


[2 

2 k [/'(£)] 2 > 


(i3) 

This relation means that the energy-momentum tensor 
in the Einstein frame is not covariantly conserved, un¬ 
less f(R) = R or T = 0 (29]. We can write the field 
equation m as 


AW - \R 9 w = «(-C + T A)> ( 14 ) 


Eqs. m and ra give the time evolution of R, 


£ V6c(Rf'-2f)y/Rf'-3f 

2 f 2 + Rf'f" - 6ff" 


(19) 


and we wrote / = f(R) for brevity. The equation of 
continuity m can be written as 


e m + 3 He m — Q, (20) 


where 


„_ kf"(R)[Rf(R)-2f(R)] 

Q 2 K[f'{R)Y 

The quantity Q describes the interaction between mat¬ 
ter and dark energy j33j, and vanishes for the general- 
relativistic case f(R) = R. 

In the comoving frame of reference, Eq. m reads 


A 

—-h 

K 


i - f\R) 
f'{R) £m ’ 


A 


PA = -• 

K 


( 22 ) 

(23) 


The definitions in Eq. were chosen so that O m -f-OA — 
e m /e c + ca/cc = 1 , where e c = 3 H 2 /{kc 2 ) is the critical 
energy density, 


where T™ = This defines the dark energy- 

momentum tensor, 




A (R) 

—— 9nv 


1 ~ f'{R) 

f'(R ) 


(15) 


Rf(R) - 3 f(R) 

Mrm 2 


(24) 


We also find the equation of continuity for dark energy, 


From Eq. lfT~Tll it follows that matter and dark en¬ 
ergy form together a system that has a conserved four- 
momentum. Therefore, we can speak of an interaction 
between matter and dark energy [3C}- In the Jordan 
frame of f{R) gravity, the energy-momentum tensor is 
always conserved |3l[ as a direct consequence of Eq. © 
for the case of an infinitesimal translation Jl|. 

For a flat .32] Robertson-Walker universe with pres¬ 
sureless matter (dust) p m = 0 and in the comoving frame 
of reference, Eq. <P1) gives the first Friedmann equation, 


f&h7 


.(*) V(R) 


+ 3H(R)e m (R)=0, (16) 


dt [ f'(R) k 
where the Hubble parameter is given by m 


H(R) = 


lRf'(R)-3f(R) 


f'(R ) 


6 


(17) 


The matter energy density e m (R) = T(R) is obtained 
from m, 


e 


m 


Rf'(R) - 2 f(R) 


e A + 3H(e A +PA) = -Q, (25) 

which has the same form as in [33]]. Eqs. m and m 
are a realization of the fact that the system matter-dark 
energy is closed. The interaction term Q corresponds to 
the production of particles from quintessence. 

In f(R) gravity, the dark energy density is given by 


ca 


Rf'{R) - 3 f{R) - 2 f'(R)[Rf(R) - 2 f(R)} 

2 4f'(W 


Consequently, the rate of interaction T = Q/e A ^3] 
equals 


V6cf”(Rf-2f) 2 y/Rf'-3f 

~ [Rf - 3/ - 2 f'{Rf - 2/)](2/'2 + Rf'f" - 6 //") ’ 

(27) 

and the nondimensional ratio 7 = T /H takes the form 

= _ 6 ff"(Rf' - 2 f) 2 _ 

7 [Rf - 3/ - 2 f{Rf - 2/)] (2/« + Rf'f" - 6ff")' 

(28) 


Kf'{R) 


(18) 
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Finally, the ratio wa = Pa/^a is given by 


where the constant b is given by 


wa = 


Rf'-f 


Rf - 3/ - 2 f(Rf - 2 /) 


(29) 


The case wa = const was examined, for example, in (341 . 
In f(R) gravity, both 7 and wa are not constant, in gen¬ 
eral. 

We now examine the consistence of the f(R) descrip¬ 
tion of interacting dark energy with some of the phe¬ 
nomenological models studied in 0 . We begin with the 
decay law A = (3H 2 /c 2 , where (3 is a nondimensional and 
positive (since A decreases with time) constant. This con¬ 
dition follows from the dimensional relations tp ~ H~ 3 
and A ~ lpf(tpi/tn) 2 , where tp is the Hubble time, and 
lpi and tpi are respectively the Planck length and time. 
It can also be deduced from Dirac’s large number hy¬ 
pothesis 0 , and the holographic principle E35j. Using 
Eqs. (1121) and 1171) we obtain 


Rf'(R)-f(R) P 

Rf'(R)-3f(R) 3’ 


(30) 


which has the solution in the form of a power function, 
f(R) ~ jf3(i+/3)/(3+/!))_ Therefore, the relation A oc H 2 
is inconsistent with f(R) gravity models which assume 
f(R) = R+sg(R), where g(R) is some nonlinear function 
of R and e is a small quantity. Such a form of f(R) is 
required for the deviations from general relativity to be 
compatible with solar system experiments |.'ld . 

Considering A oc H 2 is equivalent 0 |, at the level of 
the Einstein equations in general-relativistic cosmology, 
to the constancy of the parameter 0 


x = 


£a 


+ £A 


(31) 


i.e. the constancy of the ratio r = e m /eA 0,10. This 
is not the case in f(R) gravity, where this parameter is 
given by 


Rf - 3/ - 2 f(Rf - 2/) 
Rf - 3/ 


(32) 


and its constancy leads to a first-order differential equa¬ 
tion for f(R) which has no power-function solutions if 
x f 0 and 1 / 1 jjuj. A more general case r oc a~ 
where £ = const, was examined in 0]. The condition 
A oc H 2 has also been used in 0] , together with the con¬ 
stancy of 7 (and r). The latter leads to a complicated, 
nonlinear second-order differential equation for the func¬ 
tion f(R). A similar equation is obtained for the decay 
law 7 oc r used in |34|,0]. We conclude that f(R) grav¬ 
ity does not favor theories which constrain Q, T, or 7 . 
Lastly, we note that a more general condition A oc H n , 
where n = const f 2 , has no power-function solutions for 
f(R) as well. 

From the dependence of the redshift on R 0) we ob¬ 
tain 


3 _ b[Rf{R) - 2/(1?)] 

[f(RW 


(33) 


inRo)? 

a o[Rof(Ro) — 2 /(i? 0 )] ’ 


(34) 


and the subscript 0 denotes the corresponding present 
value. Therefore, the decay law A oc a ~ m , where to = 
const, yields 


[f(R)]^-V[f(R)-Rf(R)) 

[Rf(R)-2f{R)]V 


= 6 = const. 


(35) 


If to = 3 then Eq. iE3 has the power-function solution, 
f(R) ~ Ri 1 - 25 )/^- 6 )_ The positivity of £a and A requires 
5 < 0. For m = 2 0 , 0 ], to = 4 0 , 0 , 0 , or other 
values of to 0, Eq. ll35l) is quite complicated. 

We now express the condition A oc q [Tsf. where q is 
the deceleration parameter, in f(R) gravity. Using the 
formula for q 0 , 


we find 


q(R) = 


2 Rf(R) - 3 f(R) 

Rf{R ) - 3 f(R) ’ 


(Rf-f)(Rf-3f) 
f 2 (2Rf - 3/) 


= const. 


(36) 


(37) 


The only power function obeying this condition is f(R) = 
R 2 , which is excluded by Eq. CEB- Concluding, there is 
no simple correspondence between f(R) gravity and the 
phenomenological, power-law models examined in U5 
which relate the cosmological constant to any (one) cos¬ 
mological quantity from the set (H, q , a). In order to 
find the function f(R) which corresponds to such an ex¬ 
pression for A, we need to solve a complicated, nonlinear 
differential equation. 

However, using Eqs. CE3, 03, and (p| , we can com¬ 
pose the exact cosmological term m from any two cos¬ 
mological paremeters from the set ( H , q, a ), which gives 
three combinations. We find the following relations: 


H 2 


-r(i -2 9), 

(38) 

3 H 2 a~ 3 

(39) 

c 2 b ’ 

1 — 2 q a~ 3 

2(1 + q) b 

(40) 


All three are satisfied in f(R) gravity, but we only need 
one to compare f(R) gravity with a particular phe¬ 
nomenological model. The first equation generalizes the 
law A oc H 2 in the sense that the proportionality con¬ 
stant depends on the deceleration parameter, (3 = 1 — 2 q. 
Eq. 131 has the same form for an arbitrary function 
f(R), and is an important prediction of the Palatini f(R) 
gravity |0. The two other relations involve the constant 
b which depends on the function f(R) and the present 
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value i? 0 , i.e. they are model dependent. Eq. is 
similar to the form assumed in 1 1311 . except we have a~ 3 
instead of a" 2 . We note that if A depends on q then the 
relation q = 1/2, which holds for the matter era, yields 
A = 0. 


In general-relativistic cosmology, some phenomenologi¬ 
cal laws are dynamically equivalent to each other 4(1 l4p. 
For example, the relations A oc H 2 , A oc qH 2 [43, and 
A oc e m ;46| give the same cosmological evolution ITif . 
Since the decay law A oc H 2 appears to be inconsistent 
with current constraints on the dark energy equation of 
state 0,E3i so do the other equivalent models. In f(R) 
gravity, the condition A = qiqH 2 /c 2 , where qi = const, 
becomes 


Rf'(R)-f(R) _ ^ 
2Rf'(R)-3f(R) 3' 


appears (with Palatini variation and in the ECF) to be 
compatible with cosmological data, although more obser¬ 
vations are necessary to put stronger constraints on the 
analyzed parameters 1411. The current value of R equals 
Ro = (—1.05 ± 0.01)ck |2!j, which gives the present value 
of 7, 

70 = 0.015 ± 0.005. (44) 

Such a small value indicates that the interaction between 
matter and dark energy is weak, as compared to the direct 
interaction between dark energy and spacetime. A large 
error arises from the sensitivity of the term Rf(R) — 
2 f(R) around Rq. 

Similarly, we find the value of 7 at the deceleration- 
to-acceleration transition which occurred at Rt = 
-sJbJZa 0 : 


Eq. 1411) has the solution in the form of a power func¬ 
tion, f(R) ~ /j 3 ( 1 +/d/( 3 + 2 ^). This solution is equivalent 
to that of Eq. EB if we associate n = 2/3/(l — /3), in 
agreement with If .1 0- On the other hand, the condi¬ 
tion A oc e m leads in f(R) gravity to 


Rf'-f 

f'(Rf - 2/) 


= const, 


(42) 


which has no power-function solutions and is not equiv¬ 
alent to the law A oc H 2 . 

Finally, we should mention that the decay law (El pre¬ 
dicted by f(R) gravity resembles the form used in |47f . 
A oc R = — 6U 2 (1 — q)/c 2 . This form is also equiva¬ 
lent, in gene ral-relativistic cosmology, to the condition 
A oc II 2 Thill ■ R- e f- Si showed that almost all the cur¬ 
rent phenomenological models of decaying vacuum (dark 
energy) can be unified in the form of the modified mat¬ 
ter expansion rate, e m oc a~ 3+£ , where e = const [56|. 
The relation A a H 2 requires e > 1 to get a currently 
accelerating universe. For example, the results of l4ll| 
give £0 = 1-85 E?j. However, this inequality leads also 
to an accelerating expansion of the matter dominated 
universe. Overall, the type la supernovae, cosmic mi¬ 
crowave background, and large-scale structure observa¬ 
tions give the constraint e < 0.1, which rules out the law 
A x II' 2 00. 


IV. THE R - 1/7?. GRAVITY 

The consistence of f(R) models with cosmological data 
has been studied for both metric and Palatini variational 
formalisms 01 ■ We examine a particular case 

f{R) = R ~^ {43) 

where a is a constant, which is the simplest way to in¬ 
troduce current cosmic acceleration in f(R) gravity [f|. 
This model is referred to as the R — 1/R gravity, and 


7. = A- (45) 

The interaction rate drops significantly (~ 20 times) from 
the time of the transition to the present, which does not 
support the assumption of the constancy of 7 0, El • 
As the universe approaches asymptotically a de Sitter 
expansion, R tends to —a and the interaction rate 7 be¬ 
tween matter and dark energy decreases to zero. 

In the case of the R — 1/R gravity, Eq. (l29l) for the 
present time gives 

to A ,o =-1.10 ±0.02. (46) 

This result is consistent with the observed wa,o = 
— 1.02ig^g 0|. We also find the value of w\ at the 
deceleration-to-acceleration transition, 

WA,t = - (47) 

As the universe approaches a de Sitter phase, iva tends 
to the value —1. Such a value corresponds to dark en¬ 
ergy which does not interact with other forms of matter. 
Similarly, Eq. |32|) reads 


x 0 = 0.69l°;° 3 , (48) 

*t = \- (49) 

The observational constraint on the ratio x in the matter 
era is x < 3 x 10 -5 0]. The deceleration-to-acceleration 
transition occurred when the universe was, clearly, dom¬ 
inated by dark energy. As the universe approaches a de 
Sitter expansion, x tends to 1. 

From Eq. E3), we obtain the present value for the 
cosmological term A. Substituting the observed values 
H 0 = 71 ±4km s _1 Mpc _1 0 and q 0 = -0.74i0.18 0, 
we find 


A 0 = (l-46t°;") x 10~ 52 m -2 . 


(50) 
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If we use the gn predicted by the R — l/R gravity, an = 
-0.671°;°® 0, we arrive at 

A 0 = (1.381^2°) x 10“ 52 m“ 2 . (51) 

The cosmological constant at the deceleration-to- 
acceleration transition is given by A t = H//c 2 . Its 
value depends on the form of the function f(R). For 
the case (I43[) . we use the expressions for H tl a, and go 
found in (29j to obtain 

^=0.96 ±0.08. (52) 

A t 

We see that A has not changed much since the time when 
q changed the sign. A significant decrease of the cosmo¬ 
logical term must have ha ppe ned earlier. 

Using fi TOi o = 0.29l°'o3 [32l| then we obtain 

— = 27.211°. (53) 

7o 

This value differs significantly from that found in j44j, 
r0/70 = 0.54. For the deceleration-to-acceleration transi¬ 
tion (which is in the range of redshifts examined by §52}), 
we have fl m ,t = 4/5 This gives 

- = 15, (54) 

7 1 

which is again much larger than ~ 1. Therefore, the 
condition relating the dark energy density to the Hubble 
parameter together with the assumption of the constancy 
of the interaction rate are incompatible with the R — 1 /R 
gravity. This conclusion agrees with the results of [43 
which disfavor the decay law A oc H 2 , and the results 
of m which support f / 0. 

Eq. CD can be written as 

£m + nHe m = 0, (55) 

where 

n = 3-7(H“ 1 -l). (56) 

The R — l/R gravity predicts for the present time 

n 0 = 2.96 ±0.01. (57) 

For the moment of the deceleration-to-acceleration tran¬ 
sition, the R—l/R gravity gives 

n t = 2.93. (58) 

The departure of n from 3 means that matter alone is 
not conserved. The above values of n are smaller than 
3 (unlike 36/11 found in §21), indicating that matter is 
being produced from dark energy. It has been shown [29i| 
that the largest deviation from the standard nonrela- 
tivistic matter scaling occurs around the deceleration-to- 
acceleration transition, where ne m ~ a. The result (|55l) 
is a numerical estimation of this deviation. Since n is 


very close to 3 at this transition, and the difference n — 3 
is negligible in the early universe |29} . we may say that 
n ~ 3 for the entire matter era. Therefore, the devia¬ 
tion of the growth of the cosmic scale factor in this era 
from the standard law a(t) ~ t 2 / 3 is very small which 
is consistent with WMAP cosmological data §5(3 ■ The 
values n = 2.93 and n = 2.96 correspond respectively to 
e = 0.07 and e = 0.04, in agreement with the observa¬ 
tional constraint e < 0.1 J43j. Our results show that the 
R — 1 /R gravity in the Palatini variational formalism is 
a viable theory of gravitation that explains current cos¬ 
mic acceleration. On the other hand, the compatibility 
of metric f(R) gravity models with astronomical obser¬ 
vations is still an open problem (52] ■ 

Finally, we compare the R— 1/R gravity with the gen¬ 
eral observational constraint on the matter-dark energy 
interaction, found in l39l |. This constraint is given by 


r 

v H- 

H(l + r) 


< 7 < — 


3wac 
1 ± r 


(59) 


However, it is easier to use the conditions from which 
Eq. jSHt was derived. The first one, go < 0 (accel¬ 
eration), is satisfied by the model 14:511 which predicts 
go = — 0.67to'o3 2U . The second one represents our 
expectation that the ratio of the matter density to the 
dark energy density decreases with the evolution of the 
universe, r < 0. In f(R) gravity, this condition yields 


ft d [ Rf - 3/ 
dR[f'(Rf'-2f)_ 


> 0 , 


(60) 


where R is given by d. For the case m, the quan¬ 
tity R is positive in the range of acceleration [R t ,Roo] = 
[— \/5/3a, —a}. Therefore, we obtain 


_d_ 

dR 


p 2 a 
R --3R 


(1 


3R 2 




> 0 , 


(61) 


which holds so long as R < —a/3. The last inequality 
is satisfied for both the matter and dark energy eras of 
the universe expansion 001 , and so is the condition 
r < 0. 


V. SUMMARY 

f(R) gravity provides a relativistically covariant expla¬ 
nation for a cosmological constant that varies with time 
and interacts with matter, basing on the principle of least 
action. We analyzed the interaction between matter and 
dark energy in f(R) gravity formulated in the Einstein 
conformal frame. We used the Palatini variational prin¬ 
ciple to obtain the field equations and apply them to a 
flat, homogeneous, and isotropic universe filled with dust. 
We found how the simplest phenomenological models of a 
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variable cosmological constant are related to f(R) grav¬ 
ity. Particularly, f(R) gravity predicts, for a flat universe 
without pressure, a simple relation Ac 2 = H 2 { 1 — 2 q). 

For the particular case f(R) = R — a 2 /(3R), we found 
that the interaction rate changes significantly between 
the moment of the deceleration-to-acceleration transi¬ 
tion and now. During the same period, the cosmolog¬ 
ical constant does not change much, indicating that its 
significant decrease must have happened earlier. The pre¬ 
dicted value of the current interaction rate is on the or¬ 
der of 1% of the present value of the Hubble parame¬ 
ter, which means that this interaction is relatively weak. 
Consequently, the energy density scaling only slightly 


(by ~ 1%) deviates from the standard scaling for non- 
relativistic matter, which is consistent with cosmological 
data. Therefore, f(R) gravity in the Palatini formalism 
appears as a viable theory that explains current cosmic 
acceleration. 

All our predictions for nondimensional cosmological 
quantities in the f(R) = R — a 2 /(3 R) gravity are inde¬ 
pendent of the value of the only parameter in this model, 
a. This is not true for more complicated cases where 
a function f(R) contains two or more parameters from 
which one can compose one or more nondimensional com¬ 
binations. 
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